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Abstract 

We report improved measurements of the partial rate asymmetry (Acp) in B -» M decays 
with 140 fb _1 of data collected with the Belle detector at the KEKB e + e~ collider. Here 
h stands for a charged or neutral pion or kaon and in total five decay modes are included: 
K^ir*, K%tt*, K*tt°, vr^vr and Kgir . The flavor of the last decay mode is determined from the 
accompanying B meson. Using a data sample 4.7 times larger than that of our previous measure- 
ment, we find AcpiK^) -0.088 ±0.035 ±0.013, 2.4a from zero. Results for other decay modes 
are also presented. 

PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw, 14.40.Nd 
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In the Standard Model (SM) CP violation arises via the interference of at least two 
diagrams with comparable amplitudes but different CP conserving and violating phases. 
Two-body charmless hadronic B decays may thus be used to access CP violation because of 
the possibility of large interference between tree and penguin diagrams. The angle 0i of the 
unitarity triangle has been measured with good accuracy using b — > ccs transitions 0, HI • 
The focus of the B factory experiments has now shifted to measurement of the other two 
angles, 02 and 03, and to the search for direct CP violation. 

Recent theoretical work suggests that 02 and 03 could be constrained from measure- 
ments of branching fractions and partial rate asymmetries in B — > hh decays. However, 
different theoretical approaches such as perturbative QCD (PQCD) [3] and QCD factor- 
ization (QCDF) @, i, Q lead to different predictions for the partial rate CP asymmetry, 
defined as 

N(B -> /) ~ N(B -> f) 
CP N(B —»•/) + N(B — > /) ' {) 

where N(B — > /) is the yield for the B — > hh decay and N(B — > f) denotes that of the 
charge-conjugate mode. For instance, QCDF predicts a positive asymmetry of less than 10% 
in B — > while PQCD favors a negative value ranging from —10% to —30%. Although 
there are large uncertainties related to hadronic effects, more precise measurements of partial 
rate asymmetries would discriminate between the theoretical approaches. Furthermore, if a 
sizable asymmetry is observed for a mode where only one diagram contributes in the SM, 
this may be an indication of new physics. 

In the search reported in this paper, five decay modes are considered: 
K~^7r^~ , KgiT T , K t tt°, 7c t tt° and K^tt . (In the following, charge-conjugate modes are implied 
in the following unless otherwise stated.) The data sample used in this study corresponds to 
140 fb _1 (152 x 10 6 BB pairs), which is 1.8 times larger than that used in the previous study 
for the Kgii~ mode Q, and is 4.7 times larger for the other modes The observation of 
large CP violation in B° — ► 7r + 7r~ at Belle has been recently reported elsewhere |9(. 

The data were collected with the Belle detector at the KEKB e + e~ collider. KEKB is an 
asymmetric collider (3.5 on 8 GeV) that operates at the T(4S) resonance (y/s = 10.58 GeV). 
The Belle detector is a large-solid-angle magnetic spectrometer consisting of a three-layer 
silicon vertex detector, a 50-layer central drift chamber (CDC), a system of aerogel threshold 
Cherenkov counter (ACC), time-of-flight scintillation counters (TOF), and an array of Csl 
crystals (ECL) located inside a superconducting solenoid that provides 1.5 T magnetic field. 
An iron flux return located outside is instrumented as a K° L and muon identifier. A detailed 



description of the Belle detector can be found in Ref. [R 

The B candidate selection is the same as described in Ref. jllj. Charged tracks are 
required to originate from the interaction point (IP). Charged kaons and pions are identified 
using dE/dx information and Cherenkov light yields in the ACC. The dE/dx and ACC 
information are combined to form a K-tt likelihood ratio, TZ(Ktt) = Ck/{£>k + Ar)> where 
Ck/tt is the likelihood of kaon/pion. Charged tracks with TZ(Ktt) > 0.6 are regarded as 
kaons and 1Z(Ktt) < 0.4 as pions. The K/n identification efficiencies and misidentification 
rates are determined from a sample of kinematically identified D* + — > D°ti + ,D — > K~tt + 
decays. For positively charged tracks, the identification efficiencies and fake rates for K(tt) 
are 83.7±0.2% (91.3±0.2%) and 10.7±0.2% (5.1±0.1%) and for negatively charged tracks, 
they are 84.7 ±0.2% (90.5 ±0.2%) and 10.0 ±0.2% (5.7 ±0.1%), respectively. Furthermore, 
charged tracks that are positively identified as electrons are rejected. 
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Candidate 7r° mesons are reconstructed by combining two photons with invariant mass 
between 115 MeV/c 2 and 152 MeV/c 2 , which corresponds to ±2.5 standard deviations. Each 
photon is required to have a minimum energy of 50 MeV in the barrel region (32° < 6 1 < 
129°) or 100 MeV in the end-cap region (17° < 6 1 < 32° or 129° < 6 1 < 150°), where # 7 
denotes the polar angle of the photon with respect to the beam line. To further reduce 
the combinatorial background, 7r° candidates with small decay angles (cos#* > 0.95) are 
rejected, where 9* is the angle between a 7r° boost direction from the laboratory frame 
and its 7 daughters in the 7r° rest frame. Candidate K s mesons are reconstructed from 
pairs of oppositely charged tracks with invariant mass in the range 480 MeV/c 2 < M W7r < 
516 MeV/c 2 . Each candidate must have a displaced vertex and a flight direction consistent 
with a K s originating from the IP. 

Two variables are used to identify B candidates: the beam-constrained mass, M^ c = 



i?beam — p* B , and the energy difference, AE = E* B — E^ e3um , where E^ eaum is the beam 
energy and E* B and p* B are the reconstructed energy and momentum of the B candidates 
in the center-of-mass (CM) frame. Events with M^ c > 5.20 GeV/c 2 and —0.3 GeV < 
AE < 0.3 GeV are selected for the final analysis for modes with 7r°'s in the final state 
and —0.3 GeV < AE < 0.5 GeV for the modes without. For the K s tt° mode, events with 
AE < —0.1 GeV are excluded to remove the background from charmless B decays. 

The dominant background is from e + e _ — > qq (q = u, d, s, c) continuum events. To 
distinguish the signal from the jet-like continuum background, event topology variables and 
B flavor tagg ing information are employed. We combine a set of modified Fox- Wolfram 
moments 12j| into a Fisher discriminant. The probability density functions (PDF) for this 
discriminant, and the cosine of the angle between the B flight direction and the z axis, are 
obtained using signal and continuum Monte Carlo (MC) events. These two variables are 
then combined to form a likelihood ratio 1Z = C S /(C S + £ q q), where C s {qq) is the product 
of signal (qq) probability densities. Additional background discrimination is provided by 
B flavor tagging. The Belle standard flavor tagging algorithm [l3[ gives two outputs: a 
discrete variable indicating the flavour of the tagging B and the MC-determined dilution 
factor r, which ranges from zero for no flavor information to unity for unambiguous flavor 
assignment. An event that contains a lepton (r close to unity) is more likely to be a BB 
event so a looser 7Z requirement can be applied. We combine 1Z and r to form a multi- 
dimensional likelihood ratio (MDLR), defined as Cf DhR /(Cf DLR + £^- DLR ), where 
denotes the likelihood of the signal (qq) in the TZ-r space. The continuum background is 
reduced by applying a selection requirement on the MDLR. The requirement for each mode 



is optimized according to the figure of merit defined as N^ xp / \J Ns Xp + N^ p , where N^ xp 
denotes the expected signal yields based on our previous branching fraction measurements 
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and Nqq P denotes the expected qq yields from sideband data (M bc < 5.26 GeV/c 2 ). A 
typical requirement suppresses 94-98% of the continuum background while retaining 50-65% 
of the signal. 

Backgrounds from T(4S) — > BB events are investigated using a large MC sample. After 
the MDLR requirement, we find a small charmless three-body background at low AE, and 
reflections from other B — > hh decays due to K-n misidentification. 

The signal yields are extracted by applying unbinned two dimensional maximum likeli- 
hood (ML) fits to the (M bc and AE) distributions of the B and B samples. The likelihood 
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for each mode is defined as 

L = e —^- xniEW and (2) 

Vj = ^[1 - g< • Acpj\Pj{Mhcu A^i), (3) 

where A" is the total number of events, i is the identifier of the z-th event, P(M^ C , AE) 
is the two-dimensional PDF of M^ c and AE, q indicates the B meson flavor, B(q = +1) 
or B(q = —1), Nj is the number of events for the category j, which corresponds to either 
signal, qq continuum, a reflection due to K-tt misidentification, or background from other 
charmless three-body B decays. 

Unlike the other four hh decay modes, the flavor of the B meson in the K^ir channel is 
not self-tagged and must be determined from the accompanying B meson. To account for the 
effect of BB mixing and imperfect tagging, the term Aqp for the signal in Eq. 3 has to be 
replaced by *4 CP (l-2xd)(l-2w z ), where \d = 0.181±0.004 is the time-integrated mixing 



parameter and and wi is the wrong-tag fraction. The K^tc sample is divided into six r-bins, 
and the r-dependent wrong-tag fractions, wi (I = 1,...,6), are determined using a high 
statistics sample of self-tagged B° — > D^" n + , D*~ p + and D*~t r v events and their charge 
conjugates [15]. No K-ir misidentification reflection or other backgrounds from charmless 
three-body B decays are included in this mode. 

The yields and asymmetries for the signal and backgrounds are allowed to float in all 
modes except for Kgir , in which the background asymmetry is set to zero. Since the 
K + 7t° and tt + tt° reflections are difficult to distinguish with AE and M^ c , we fit these two 
modes simultaneously with a fixed reflection-to-signal ratio based on the K-n identification 
efficiencies and fake rates. All the signal PDFs (P(M^ C , AE)) are obtained using MC sim- 
ulation. No strong correlations between M^ c and AE are found for the K° s ir~ and K~ir + 
signals. Therefore, their PDFs are modeled by products of single Gaussians for My, c and 
double Gaussians for AE. For the modes with 7r°'s in the final state, there are correla- 
tions between M^ c and AE in the tails of the signals; hence, their PDFs are described by 
smoothed two-dimensional histograms. Discrepancies between the signal peak positions in 
data and MC are calibrated using B + — > D°n + decays, where the D° — > K + 7!-~n° sub-decay 
is used for the modes with a 7r° meson while D° — > K + tt~ is used for the other modes. The 
MC-predicted AE resolutions are verified using the invariant mass distributions of high mo- 
mentum D mesons. The decay mode D° — > K + it~ is used for B° — > K + tt~, D + —* Kg7r + for 
B + — > Kg7r + and D° — > K + it~tt for the modes with a 7r° in the final state. The parameters 
that describe the shape of the PDFs are fixed in all of the fits. 

The background PDFs for continuum events are modeled by an ARGUS function jl6| for 
the Mb c distribution and a 1st or 2nd order polynomial for the AE distribution, assuming 
no correlations between the two. A large MC sample is used to investigate the background 
from charmless B decays and a smoothed two-dimensional histogram is taken as the PDF. 

Table E gives the signal yields and Aqp values for each mode. The asymmetries for the 
background components are consistent with zero. The errors on Aqp are determined from 
the ML fit and are consistent with MC expectations. Projections of the fits are shown 
in Fig. ^ The systematic errors from fitting are estimated by checking the deviations 
of the Acp after varying each parameter of the signal PDFs by one standard deviation. 
The uncertainty in modeling the three-body background is studied by excluding the low 
AE region (< —0.12 GeV) and repeating the fit. For the K + ti and 7r + 7r° modes, the 
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FIG. 1: Mb c (left) and AE (right) distributions for B (left two columns) and B (right two columns) 
candidates reconstructed in the K^ir^, Kgir^, K^tt , tt^tt and K^tt modes (from top to bottom). 
The histograms represent the data, while the curves represent the various components from the 
fit: signal (dashed), continuum (dotted), three-body B decays (dash-dotted), background from 
mis- identification (hatched) , and sum of all components (solid) . 



uncertainty on the reflection-to-signal ratios obtained from D* control samples are included 
in the systematic error, with contributions of ±0.006 and ±0.009 for the K + tt° and 7r + 7r° 
modes, respectively. The uncertainties in the mixing parameter Xd and wrong tag fraction 
w r are included in the systematic error for the K^n mode. At each step, the deviation in 
Acp is added in quadrature to provide the systematic errors due to fitting, which are io!o20 
for the KgiT mode and less than 0.01 for all other modes. A possible detector bias in Acp 
is investigated using the B + — > D°tt + samples, assuming null charge asymmetries in the 
production. In this control sample, the same continuum suppression and fitting procedures 
are applied. To mimic the B° — > K^tt analysis, the flavor of the B meson is determined 
using flavor tagging information from the associated B. Results of these null asymmetry 
checks are shown in Table ITT1 The final systematic errors are then obtained by quadratically 
summing the errors from the null asymmetry tests and the fitting systematics. 

In this study, the partial rate asymmetry Acp{K~tc + ) is found to be —0.088 ± 0.035 ± 
0.013, which is 2.4<r from zero. The corresponding 90% confidence level (C.L.) interval is 
—0.15 < Acp(K~n + ) < —0.03. Our central value is similar to that reported by BaBar, 
Acp(K~n + ) = — 0.107±0. 041 ±0.013 |17| . indicating that the partial rate asymmetry may 
be negative. Theoretical predictions from different approaches suggest that Acp(K 
and Acp{K + 7t°) should have the same sign. The uncertainty in our result for Acp{K + tt°), 
is large enough for it to be consistent with this expectation. We set a 90% C.L. interval of 
—0.04 < Acp(K + 7c°) < 0.16. Since no evidence of direct CP violation is observed in the 
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TABLE I: The fitted signal yields and Aqp results for individual modes. 



Mode Signal Yield A C p 



K T 7r ± 1029.8 ± 35.3 

K^ir? 429.8 ± 23.0 

K^tt 516.5 ±29.4 

tt^tt 235.6 ± 23.4 

K%ir 96.6 ± 11.8 



-0.088 ± 0.035 ± 0.013 
+0.05 ±0.05 ±0.01 
+0.06 ± 0.06 ± 0.02 
-0.00 ±0.10 ±0.02 
+0.16 ±0.29 ±0.05 



TABLE II: Results of null asymmetry cross-checks. 

Mode B -► £>(-► Ktt)tt B -► D{-^ Kirir°)ir B ^ Dir tagged 
N sig 6290.2 + 79.8 7982.1 ± 101.8 8726.1 ±86.4 
A C p 0.002 ±0.013 0.012 ±0.012 -0.027 ± 0.029 



decays B~ -> K° s n~ and B + -> tt+tt , we set 90% C.L. intervals: -0.04 < .AcpC^") < 
0.13, -0.17 < Acp(n + n°) < 0.16 and -0.33 < A C p{K° s 7r ) < 0.64. 
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